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The  ability  to  successfully  perform  many  aspects  of  a  satellite  mission  is  directly  im¬ 
pacted  by  the  ability  to  precisely  determine  and  predict  the  satellite’s  orbit  through  high 
precision  orbit  determination.  The  orbit  determination  process  relies  on  numerical  pro¬ 
cedures,  satellite  measurements,  and  force  modeling  to  generate  the  orbit  solution  and 
prediction.  As  a  foundation  for  detailed  force  modeling,  gravity  forces  due  to  the  distri¬ 
bution  of  the  Earth’s  mass  must  be  accurately  modeled.  In  prior  research,  a  number  of 
geopotential  models  and  ocean  tide  models  have  been  developed  for  use  in  orbit  determi¬ 
nation.  This  paper  examines  current  and  historically  recommended  geopotential  and  ocean 
tide  models  using  the  Naval  Research  Laboratory’s  Orbit  Covariance  Estimation  and  ANal- 
ysis  tool.  Geodetic  satellites  with  high  precision  satellite  laser  ranging  measurements  are 
used  as  test  cases  to  evaluate  the  predictive  capabilities  of  the  geopotential  and  ocean  tide 
models.  Orbit  fit  and  prediction  consistency  metrics  are  generated  for  multiple  geopotential 
and  ocean  tide  model  combinations.  Results  show  that  use  of  the  EGM2008  geopotential 
model  and  the  GOT4.8  or  FES2004  ocean  tide  models  generally  result  in  predictive  orbit 
solutions  that  more  closely  follow  the  definitive  orbit  solution.  However,  these  results  vary 
for  different  satellite  orbits  and  time  past  the  initial  fit  span. 

I.  Introduction 

The  ability  to  perform  a  range  of  satellite  missions  is  directly  impacted  by  the  ability  to  precisely  deter¬ 
mine  and  predict  a  satellite’s  orbit.  Detailed  force  modeling,  including  gravity  forces  due  to  the  distribution 
of  the  Earth’s  mass,  must  be  included  in  the  orbit  determination  process  to  accurately  determine  and  predict 
the  satellite  orbit.  A  number  of  geopotential  models  and  ocean  tide  models  have  been  developed  and  may 
be  used  to  model  gravity  forces  to  varying  degrees  of  success.  The  Naval  Research  Laboratory  (NRL)  Orbit 
Covariance  Estimation  and  ANalysis  (OCEAN)  tool  is  used  to  evaluate  the  suitability  of  these  models  to 
the  application  of  orbit  determination. 

OCEAN  is  a  highly  configurable,  database  driven  software  tool  that  enables  precision  orbit  determination 
for  a  range  of  satellite  missions.  OCEAN  allows  users  to  simulate  data,  propagate  a  spacecraft  state,  or  solve 
for  an  orbit  using  a  Kalman  Filter-Smoother  (KFS)  or  Weighted  Least  Squares  Orbit  Determination  (WLS- 
OD)  process.  Early  history  of  OCEAN  is  given  in  Reference  1,  while  references  2,  3  and  4  discuss  further 
developments.  More  recently  OCEAN  has  been  used  to  calculate  orbits  to  support  operations  for  the  NRL 
UPPERSTAGE  and  TACSAT-4  satellite  missions. 

Previous  work  evaluates  the  suitability  for  older  geopotential  and  ocean  tide  models  in  precision  orbit 
determination5.6  In  this  paper,  the  currently  recommended  geopotential  and  ocean  tide  models  are  evaluated 
along  with  historically  recommended  models  to  better  understand  expected  performance. 

II.  Orbit  Determination  Methodology 

The  International  Laser  Ranging  Service  (ILRS)  provides  global  satellite  laser  ranging  data  in  support  of 
geodetic  research  activities.7  The  ILRS  catalogs  Satellite  Laser  Ranging  (SLR)  data  to  a  number  of  geodetic 
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satellites.  These  satellites  are  typically  designed  to  facilitate  the  study  of  Earth’s  gravity  and  contain  laser 
retro-reflectors  to  facilitate  laser  ranging  data  collection.  High  precision  laser  ranging  measurements  can 
then  be  used  to  perform  precision  orbit  determination.  ILRS  guidelines  call  for  a  precision  of  the  normal 
point  laser  range  measurement  to  LAGEOS-1  of  under  one  centimeter.8 

OCEAN  is  used  with  SLR  data  to  calculate  precision  orbits  for  several  geodetic  satellites  including 
LAGEOS-1,  LAGEOS-2,  STARLETTE  and  STELLA.  By  comparing  predicted  orbits  to  fitted  orbits  for 
each  satellite  using  a  combination  of  geopotential  and  ocean  tide  models,  the  predictive  consistency  of  the 
underlying  models  may  be  evaluated.  Orbits  are  determined  with  the  OCEAN  Weighted  Least  Squares 
Orbit  Determination  (WLS-OD)  methodology  using  successive  five  day  increments  of  SLR  data.  The  orbit 
solution  from  the  first  five  day  data  arc  is  propagated  forward  in  time  to  thirty  days.  The  WLS-OD  process 
is  repeated  for  successive  five  day  data  arcs.  These  orbit  solutions  are  then  compared  to  the  predicted  orbit 
from  the  first  data  arc  solution.  Thirty  days  was  chosen  as  a  comparison  time  span  to  demonstrate  the 
longer  term  variation  in  predictive  consistency.  The  degree  to  which  the  predicted  and  fitted  solutions  agree 
will  be  used  as  a  metric  to  evaluate  the  use  of  the  geopotential  and  ocean  tide  models  in  precision  orbit 
determination.  Consistency  (agreement)  between  predicted  and  fitted  orbits  using  a  particular  geopotential 
and  ocean  tide  model  combination  indicates  that  the  force  models  are  capturing  the  dynamics  of  the  real 
system.  The  prediction/fit  comparison  methodology  is  depicted  in  Figure  1. 

Figure  1:  Depiction  of  Orbit  Solution  Comparison  Methodology 
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OCEAN  is  further  used  to  fit  an  orbit  solution  to  thirty  days  of  data.  The  RMS  of  the  error  residual  is 
used  as  a  metric  to  compare  the  relative  effectiveness  of  the  geopotential  and  ocean  tide  models  to  capture 
the  long  term  variation  in  the  satellite’s  orbit.  While  data  editing  methodologies  and  raw  data  quality  often 
have  a  large  impact  on  the  RMS  of  the  error  residuals,  it  is  nonetheless  used  as  a  metric  to  compare  orbit 
solution  fit  quality  across  the  various  models. 

The  OCEAN  WLS-OD  capability  employs  extensive  spacecraft,  measurement,  and  force  modeling  to 
estimate  the  desired  spacecraft  state  and  parameters.  Force  models  include  solid  Earth  tides;  pole  tides;  lunar 
and  solar  third  body  gravitational  effects;  indirect  lunar  oblateness;  general  relativistic  effects;  atmospheric 
variability;  drag;  and  solar  radiation  pressure.  OCEAN  WLS-OD  accounts  for  various  time  systems,  including 
TAI  and  UTC  time.  The  WLS-OD  functionality  also  uses  Earth  Orientation  Parameter  (EOP)  models 
from  the  International  Earth  Rotation  and  Reference  Systems  Service  (IERS),  which  account  for  precession 
and  nutation,  Earth  rotation,  and  polar  motion.  For  this  application,  a  ninth  order  multi-step  predictor- 
corrector  algorithm  is  used  to  perform  the  integration  of  the  state  variables  and  state  transition  matrix. 
OCEAN  follows  the  standards  described  in  the  IERS  2010  Conventions.9  OCEAN  also  performs  iterative 
data  editing  to  prevent  low  quality  data  from  biasing  the  orbit  solution. 

A.  LAGEOS-1 

The  LAGEOS-1  satellite  was  launched  in  1976  by  the  National  Aeronautics  and  Space  Administration 
(NASA)  as  part  of  the  Earth  and  Ocean  Dynamics  application  program.  It  was  designed  to  provide  a  long- 
lasting  laser  target  in  a  well  defined  orbit.  The  LAGEOS-1  satellite  enabled  researchers  to  study  a  range  of 
geophysical  phenomena  with  improved  accuracy,  including  the  Earth’s  geopotential10.11  The  satellite’s  low 
ballistic  coefficient  combined  with  its  spherical  shape  minimize  the  orbital  uncertainty  due  to  drag  and  solar 
radiation  forces. 
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The  nominal  orbital  elements  are  given  in  Table  1. 


Table  1:  Nominal  Orbital  Elements  for  LAGEOS-1 


Element 

Nominal  Value 

Semi-major  Axis 

Eccentricity 

Inclination 

12,  240  km 

0.0045 

109.84° 

B.  LAGEOS-2 

The  LAGEOS-2  satellite  was  launched  in  1992  by  NASA  and  the  Agenzia  Spaziale  Italiana.  The  LAGEOS-2 
satellite  is  almost  identical  in  design  to  the  LAGEOS-1  satellite;  however,  its  orbit  was  selected  to  provide 
increased  coverage  over  seismically  active  areas  such  as  the  Mediterranean  and  California.12 
The  nominal  orbital  elements  for  LAGEOS-2  are  given  in  Table  2. 

Table  2:  Nominal  Orbital  Elements  for  LAGEOS-2 


Element 

Nominal  Value 

Semi-major  Axis 
Eccentricity 

Inclination 

11,998  km 

0.0135 

52.64° 

C.  STARLETTE 

The  STARLETTE  satellite  was  launched  in  1975  by  Centre  Nationale  d’Etudes  Saptiales  (CNES).  The 
spacecraft  was  designed  to  improve  the  geopotential  model  and  to  study  solid  Earth  tides,  ocean  tides,  and 
polar  motion.13  STARLETTE  was  also  the  first  spacecraft  to  be  entirely  covered  by  laser  corner  reflectors 
which  allow  passive  SLR  observation  capabilities.14  The  STARLETTE  orbit  is  highly  sensitive  to  temporal 
and  zonal  variations  in  the  gravity  field. 

The  nominal  orbital  elements  for  STARLETTE  are  given  in  Table  3. 

Table  3:  Nominal  Orbital  Elements  for  STARLETTE 


Element 

Nominal  Value 

Semi-major  Axis 

Eccentricity 

Inclination 

7, 190  km 

0.0206 

49.83° 

D.  STELLA 

The  STELLA  satellite  was  launched  in  1993  by  CNES,  and  is  virtually  identical  to  the  STARLETTE  satellite. 
However,  it  provides  additional  coverage  over  the  polar  regions  due  to  its  inclination.  As  with  the  LAGEOS 
satellites,  both  STELLA  and  STARLETTE  are  spherically-shaped  spacecraft  with  low  ballistic  coefficients 
to  minimize  orbital  uncertainty  caused  by  drag  and  solar  radiation  pressure  forces.15 
The  nominal  orbital  elements  for  the  STELLA  satellite  are  given  in  Table  4. 

III.  Geopotential  Models 

The  International  Center  for  Global  Earth  Models  (ICGEM),  which  is  a  part  of  the  International  Asso¬ 
ciation  of  Geodesy  (IAG),  catalogs  Earth  geopotential  models.  The  catalog  contains  over  130  models  which 
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Table  4:  Nominal  Orbital  Elements  for  STELLA 


Element 

Nominal  Value 

Semi-major  Axis 

Eccentricity 

Inclination 

7, 178  km 

0.0206 

98.6° 

have  been  developed  from  as  early  as  1966  to  the  present  day.  Several  geopotential  models  have  been  used 
historically  for  precision  orbit  determination.  The  International  Earth  Rotation  and  Reference  Systems  Ser¬ 
vice  (IERS)  serves  the  astronomical  and  geodetic  communities  by  providing  relevant  data  and  maintaining 
standards  and  best  practices.  The  IERS  releases  technical  notes  on  a  regular  basis  that  define  best  practices 
for  Earth  rotation  computation,  modeling  and  analysis.  Historically,  the  IERS  has  recommended  the  use  of 
the  Earth  Gravitational  Model  1996  (EGM96)  geopotential16  and  currently  recommends  the  use  of  the  Earth 
Gravitational  Model  2008  (EGM2008)  geopotential  model.9  The  EGM2008  model  was  developed  using  data 
from  the  Gravity  Recovery  and  Climate  Experiment  (GRACE)  as  well  as  ocean  altimetry  and  surface  gravity 
data.17 

Additional  geopotential  models  have  been  developed  using  data  obtained  from  the  GRACE  mission. 
Principal  investigators  for  GRACE  at  the  University  of  Texas  at  Austin  Center  for  Space  Research  (CSR) 
have  developed  the  GRACE  Intermediate  Field  Model,  GIF48,18  which  was  developed  using  GRACE  data 
as  well  as  ocean  altimetry  and  surface  gravity  data.18  Researchers  at  CSR  have  also  developed  the  Grace 
Gravity  Model  (GGM)  geopotential  series.  The  first  GGM  series  model  (GGM01)  was  released  in  2003. 19 
Updates  to  the  GGM01  model  were  made  in  GGM0220  and  GGM03,21  released  in  2004  and  2008,  respectively. 
The  latest  in  the  GGM  series  available  from  ICGEM  is  GGM05S.  This  model  was  developed  using  GRACE 
data  only.22 

Further,  researchers  from  the  German  Research  Center  for  Geosciences  and  the  Groupe  de  Recherche  de 
Geodesie  Spatiale  have  developed  a  series  of  gravity  models  under  the  European  Improved  Gravity  model 
of  the  Earth  using  New  methods  (EIGEN)  name.  This  model  series  makes  use  of  CHAMP,  GRACE  and 
LAGEOS  satellite  data,  as  well  as  ocean  altimetry  and  surface  gravity  data.  The  most  recent  is  a  static 
model  named  EIGEN6C3S.23 

Five  geopotential  models  are  used  in  this  evaluation:  EGM96,  EGM2008,  GIF48,  GGM05S  and  EIGEN6C3S. 

IV.  Ocean  Tide  Models 

The  geopotential  models  capture  gravity  forces  resulting  from  the  distribution  of  Earth’s  solid  mass. 
Satellite  orbits  are  additionally  affected  by  gravity  forces  resulting  from  the  distribution  of  Earth’s  ocean 
mass.  The  IERS  has  historically  recommended  the  use  of  the  Center  for  Space  Research  CSR3.0  ocean 
tide  model16,24  and  currently  recommends  the  Finite  Element  Solution  2004  (FES2004)  model.9  FES2004 
was  produced  by  Legos  and  CLS  Space  Oceanography  Division  and  distributed  by  Aviso,  with  support 
from  CNES. 25,26  Both  of  these  models  have  been  incorporated  into  OCEAN  to  maintain  adherence  to  IERS 
recommendations  for  ocean  tide  modeling. 

Other  ocean  tide  models  have  also  been  developed  from  various  satellite  data.  The  CSR4.0  model  repre¬ 
sents  and  incremental  improvement  to  the  CSR3.0  model.27  The  Global  Ocean  Tide  Model  4.8  (GOT4.8)  is 
the  latest  in  a  series  of  models  that  use  TOPEX/Poseidon  and  other  satellite  altimetry  data,28  and  improves 
over  GOT4.7  due  to  better  processing  of  the  dry  component  of  the  troposphere.  Both  of  these  models  have 
been  implemented  within  OCEAN. 

Four  ocean  tide  models  will  be  used  in  this  evaluation:  CSR3.0,  CSR4.0,  FES2004,  GOT4.8. 

V.  Testing  Results 

As  can  be  anticipated,  when  examining  the  combined  performance  for  five  geopotential  models  and  four 
ocean  tide  models  in  determining  orbits  for  four  satellites,  a  large  number  of  results  are  generated.  Presenting 
these  data  in  a  tractable  manner  presents  a  challenge.  Since  the  EGM96  and  EGM2008  geopotential  models 
are  recommended  by  the  IERS  technical  conventions,  the  results  presented  are  focused  on  these  models. 
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However,  select  results  for  other  geopotential  models  will  be  presented  as  necessary  to  draw  additional 
conclusions. 

Separate  results  are  presented  for  each  satellite  considered.  For  each  satellite,  RMS  residual  errors  are 
given  for  each  fit  span  as  a  measure  of  orbit  solution  quality.  In  addition,  RSS  errors  between  the  orbit 
prediction  and  the  current  interval  orbit  solution  are  given  to  show  the  predictive  consistency  of  the  specific 
geopotential  and  ocean  tide  model  combination.  Finally,  plots  of  average  daily  RSS  position  difference 
between  the  predicted  and  fit  orbits  are  given  for  select  cases  as  a  metric  to  evaluate  the  consistency  of  the 
underlying  geopotential  and  ocean  models  for  precision  orbit  determination. 

A.  LAGEOS-1  Results 

Table  5  shows  the  RMS  of  the  residuals  for  a  30-day  orbit  solution  given  the  specified  geopotential  model 
and  ocean  tide  model.  This  may  be  used  as  a  general  comparison  to  determine  relative  performance  between 
the  competing  models. 

Table  5:  LAGEOS-1  RMS  Error  Residuals  for  30-day  Orbit  Fit  in  Meters 


Geopotential  Model 

CSR3.0 

CSR4.0 

FES2004 

GOT4.8 

EGM96 

0.4835 

0.4840 

0.4884 

0.5457 

EGM2008 

0.5038 

0.5048 

0.5100 

0.5629 

GIF48 

0.5132 

0.5143 

0.5198 

0.5720 

GGMS05S 

0.4902 

0.4910 

0.4958 

0.5512 

EIGEN6C3S 

0.4877 

0.4884 

0.4932 

0.5494 

As  can  be  seen  in  Table  5,  the  lowest  residuals  for  the  30-day  orbit  solution  are  found  when  the  EGM96 
geopotential  and  the  CSR3.0  ocean  tide  model  are  used.  However,  these  models  do  not  produce  most 
self-consistent  orbit  prediction  when  subsequent  five  day  fits  are  examined. 


Figure  2:  LAGEOS-1  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  (Left) 
and  LAGEOS-1  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  For  GOT4.8 
Ocean  Tide  Model  and  Various  Geopotential  Models  (Right) 
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Next,  Figure  2 


(a)  shows  the  average  daily  RSS  difference  between  the  orbit  predicted  in  the  first  five 
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day  fit  compared  to  the  orbit  solution  during  the  subsequent  five  day  fits.  It  can  be  seen  in  Figure  2  (a) 
that  both  the  EGM96  and  EGM2008  geopotential  models  have  substantially  lower  error  when  used  with  the 
GOT4.8  ocean  tide  model. 

To  demonstrate  the  impact  of  the  geopotential  model,  results  are  presented  for  all  the  considered  geopo¬ 
tential  models  used  with  the  GOT4.8  ocean  tide  model  in  Figure  2  (b).  It  can  be  seen  in  Figure  2  (b)  that 
there  is  little  difference  between  the  orbit  prediction  quality  using  any  of  the  considered  geopotential  models 
for  a  fixed  ocean  tide  model. 


B.  LAGEOS-2  Results 

Table  6  shows  the  RMS  of  the  residuals  for  a  30-day  orbit  solution  given  the  specified  geopotential  model 
and  ocean  tide  model. 


Table  6:  LAGEOS-2  RMS  Error  Residuals  for  30-day  Orbit  Fit 


Geopotential  Model 

CSR3.0 

CSR4.0 

FES2004 

GOT4.8 

EGM96 

0.4937 

0.4935 

0.4939 

0.4689 

EGM2008 

0.5855 

0.5856 

0.5859 

0.5587 

GIF48 

0.6153 

0.6156 

0.6157 

0.5878 

GGMS05S 

0.5306 

0.5304 

0.5309 

0.5048 

EIGEN6C3S 

0.5087 

0.5087 

0.5089 

0.4825 

As  can  be  seen  from  Table  6,  the  30-day  orbit  solution  residuals  are  the  lowest  for  the  EGM96  geopotential 
when  used  with  the  GOT4.8  ocean  tide  model.  However,  these  models  do  not  produce  the  lowest  difference 
between  orbit  prediction  and  subsequent  five  day  fits. 


Figure  3:  LAGEOS-2  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  (Left) 
and  LAGEOS-2  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  for  FES2004 
Ocean  Tide  Model  and  Various  Geopotential  Models  (Right) 


(a) 


(b) 


*-•  EGM96_FE52004 
*-•  EGM96_GOT48 
■  EGM96CSR3.0 
*  i  EGM96  CSR4.0 
i— +  EGM20Q8_FE52004 
+  +  EGM  2  DO  8_GOT 4  8 
EG  M  2  00  8  TOPEX3 
EGM2008  TOPEX4 


a: 

5  15 


k  10 


RSS  Position  Difference  vs  Time 


Next,  Figure  3  (a)  shows  the  average  daily  RSS  position  difference  between  the  orbit  predicted  from  the 
first  five  day  fit  compared  to  the  orbit  solution  during  the  subsequent  five  day  fits.  As  can  be  seen  in  Figure 
3  (a),  there  is  no  appreciable  difference  between  the  orbit  prediction  quality  of  the  considered  geopotential  or 
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ocean  tide  models.  The  GOT4.8  ocean  tide  model  provides  the  lowest  prediction-fit  difference  for  intervals  2 
and  3  (Days  5-15  in  Figure  1).  However,  the  FES2004  ocean  tide  model  provides  slightly  better  prediction-fit 
consistency  for  prediction  intervals  4,  5  and  6.  Figure  3  (b)  shows  that  there  is  little  difference  between  the 
orbit  prediction  quality  for  the  different  geopotential  models  when  using  the  FES2004  ocean  tide  model. 

C.  STARLETTE  Results 

Table  7  shows  the  RMS  of  the  residuals  for  a  30-day  orbit  solution  given  the  specified  geopotential  model 
and  ocean  tide  model. 

Table  7:  STARLETTE  RMS  Error  Residuals  for  30-day  Orbit  Fit 


Geopotential  Model 

CSR3.0 

CSR4.0 

FES2004 

GOT4.8 

EGM96 

3.8037 

3.7848 

3.8017 

3.8214 

EGM2008 

4.2087 

4.1852 

4.2059 

4.2133 

GIF48 

4.2841 

4.2606 

4.2818 

4.2850 

GGMS05S 

4.0977 

4.0746 

4.0945 

4.1076 

EIGEN6C3S 

3.9332 

3.9119 

3.9308 

3.9472 

As  can  be  seen  in  Table  3,  the  minimum  RMS  residual  error  for  the  30-day  orbit  solution  occurs  when 
the  EGM96  geopotential  and  the  CSR4.0  ocean  tide  model  are  used.  The  RMS  errors  are  higher  in  than 
the  LAGEOS-1  and  LAGEOS-2  test  cases.  The  orbit  of  STARLETTE  was  designed  to  be  sensitive  to 
tidal  variations  as  well  as  temporal  variations  in  Earth’s  gravity.  As  a  result,  it  is  more  difficult  to  capture 
STARLETTE’s  dynamics  over  a  larger  time  span.  However,  the  orbital  dynamics  are  more  easily  captured 
for  shorter  time  spans;  and  so,  the  RMS  error  for  the  five  day  fits  are  approximately  0.5  meters.  Therefore, 
the  predicted  orbit  versus  fit  orbit  comparison  is  still  an  applicable  metric. 

Next,  Figure  4  (a)  shows  the  average  daily  RSS  between  the  orbit  predicted  from  the  first  five  day  fit 
compared  to  the  orbit  solution  during  the  subsequent  five  day  fits. 

Figure  4:  STARLETTE  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  (Left) 
and  STARLETTE  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  for  GOT4.8 
Ocean  Tide  Model  and  Various  Geopotential  Models  (Right) 


•  EGM96_FES2004 
• — •  EGM96_GOT48 

■  EGM96_C5R3.D 
*-*  EGM96_  C5R4.0 
—  EG  M  2  GO  S_FE5 2004 
♦  ♦  EG  M  2  00  8_GOT  48 
EGM2QQ8TGPEX3 

*  EGM2008  TOPE X 4 


(a) 

RSS  Position  Difference  vs  Time 


Jj  40 
£ 
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As  can  be  seen  in  Figure  4  (a),  there  is  some  variation  in  the  RSS  orbit  difference  between  predicted 
and  fitted  orbits  for  the  considered  geopotential  and  ocean  tide  models.  The  combination  of  the  EGM2008 
geopotential  with  the  GOT4.8  tide  model  provides  the  lowest  prediction-fit  orbit  difference  through  most 
of  the  prediction  intervals.  However,  at  the  final  interval  the  best  predictive  performance  is  given  by  the 
combination  of  the  EGM96  geopotential  with  the  FES2004  ocean  tide  model. 

Figure  4  (b)  shows  the  RSS  orbit  difference  between  the  predicted  orbit  and  proceeding  five  day  orbit  fits 
for  the  GOT4.8  ocean  tide  model  and  all  the  considered  geopotential  models.  As  can  be  seen,  the  differences 
between  the  various  geopotential  models  are  relatively  small.  However,  either  EGM96  or  EGM2008  produce 
orbits  with  the  lowest  RMS  difference  depending  on  the  fit  span. 

D.  STELLA  Results 

Table  8  shows  the  RMS  of  the  residuals  for  a  30-day  orbit  solution  given  the  specified  geopotential  model 
and  ocean  tide  model. 


Table  8:  STELLA  RMS  Error  Residuals  for  30-day  Orbit  Fit 


Geopotential  Model 

CSR3.0 

CSR4.0 

FES2004 

GOT4.8 

EGM96 

5.5905 

5.6136 

5.6479 

5.7992 

EGM2008 

6.1335 

6.1487 

6.1670 

6.3505 

GIF48 

6.1125 

6.1284 

6.1470 

6.3286 

GGMS05S 

6.1408 

6.1555 

6.1735 

6.3586 

EIGEN6C3S 

6.1295 

6.1458 

6.1685 

6.3509 

Figure  5:  STELLA  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  (Left)  and 
STELLA  Average  Daily  RSS  Orbit  Difference  Between  Predicted  Orbit  and  Fit  Orbit  For  FES2004  Ocean 
Tide  Model  and  Various  Geopotentials  (Right) 


(a) 


(b) 


400 


350 


I  250 


;  iso- 


RSS  Position  Difference  us  Time 


•  •  EG  M9  6_F  E  S2 0 04 
*-*  EGM96_GQT48 

E6M96_CSR3.0 
*  EGM96_CSR4,0 
►  EGM2008  FES2G04 

♦  ♦  EGM2G0B_GOT4B 
— r  EGM2008yrOPEX3 
—  EG  M  2  00  8TOPEX4 


400 


350 


300 


;  250 


a  200 


;  iso 


100 


50 


RSS  Position  Difference  V5  Time 


EGM96FES20Q4 
EGM2008_F£S2004 
EIGEN6C3^FES2004 
GGM05SFES2004 
GIF48  FES2004 


As  can  be  seen  in  Table  8,  the  lowest  RMS  error  is  given  by  the  EGM96  geopotential  and  the  CSR3.0 
ocean  tide  model.  In  addition,  the  30-day  orbit  solution  residuals  are  much  higher  than  for  either  LAGEOS-1 
or  LAGEOS-2.  This  is  largely  due  to  the  fact  that  the  STELLA  orbit  was  chosen  to  be  particularly  sensitive 
to  the  solid  Earth  and  ocean  tides  as  well  as  the  temporal  variations  of  the  geopotential.  Thus,  the  longer 
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fit  orbit  solution  is  not  able  to  sufficiently  model  the  sensitivity  to  these  time  varying  forces,  causing  larger 
residuals.  However,  it  is  seen  that  with  the  shorter  five  day  orbit  fits,  the  RMS  residuals  are  much  lower  at 
approximately  0.5  meters  RMS  or  less.  Since  the  dynamics  of  STELLA  are  more  readily  captured  by  the 
shorter  data  span,  the  comparison  of  successive  five  day  orbit  fits  is  still  useful.  As  in  previous  cases,  the  set 
of  models  that  minimizes  the  predict-fit  orbit  difference  for  the  STELLA  orbit  is  not  the  same  as  the  models 
that  minimize  the  30-day  residuals. 

Figure  5  (a)  shows  the  average  daily  RSS  between  the  orbit  predicted  from  the  first  five  day  fit  compared 
to  the  orbit  solution  during  the  subsequent  five  day  fits.  As  can  be  seen  in  Figure  5  (a),  the  EGM2008 
geopotential  model  performs  much  better  than  the  EGM96  geopotential  model  regardless  of  the  ocean 
tide  model  used.  There  is  only  a  small  variation  in  orbit  prediction  consistency  caused  by  the  ocean  tide 
model.  However,  over  most  prediction  intervals,  the  combination  of  EGM2008  and  FES2004  gives  the  lowest 
prediction-fit  orbit  difference. 

Figure  5  (b)  shows  the  average  daily  RSS  error  between  the  predicted  orbit  from  the  first  five  day  orbit 
solution  and  the  proceeding  five  day  orbit  solutions  using  the  FES2004  ocean  tide  model  and  all  of  the 
considered  geopotential  models.  As  can  be  seen  in  Figure  5  (b),  all  geopotential  models  provide  an  orbit 
solution  with  similar  quality  prediction,  except  for  the  EGM96  geopotential  model,  which  preforms  much 
worse  than  the  other  models.  The  EGM2008  geopotential  model  shows  slightly  better  predictive  consistency 
than  GIF48,  GGM05S,  and  EIGEN6C3  through  most  intervals.  However,  in  the  final  interval,  EIGEN6C3S 
provides  better  performance  than  EGM2008. 

VI.  Conclusions 

There  are  a  number  of  results  between  the  four  satellites,  four  ocean  tide  models  and  five  geopotential 
models  considered.  From  these  results,  there  are  some  conclusions  that  may  be  drawn. 

In  the  case  of  LAGEOS-1,  it  is  seen  that  the  use  of  the  GOT4.8  ocean  tide  model  produces  orbit 
predictions  20  to  30  days  ahead  in  time  which  are  more  consistent  with  SLR  data  than  orbits  produced  using 
other  tide  models.  The  variation  in  predictive  consistency  due  to  the  geopotential  model  is  quite  small  when 
the  GOT4.8  ocean  tide  model  is  used. 

In  the  case  of  LAGEOS-2,  there  is  very  little  variation  in  the  predictive  performance  of  the  orbit  solutions 
regardless  of  the  models  used. 

In  the  case  of  STARLETTE,  the  combination  of  the  EGM2008  geopotential  and  the  GOT4.8  has  the  best 
predictive  performance  over  the  fourth  and  fifth  intervals,  where  the  orbit  difference  is  the  most  pronounced. 

In  the  case  of  STELLA,  it  can  be  seen  that  the  EGM96  geopotential  has  the  poorest  predictive  perfor¬ 
mance.  There  is  little  variation  when  any  other  geopotential  model  is  used  with  any  ocean  tide  model. 

These  results  show  that,  for  these  test  cases,  the  GOT4.8  and  FES2004  ocean  tide  models  along  with  the 
EGM2008  geopotential  generally  result  in  predictive  orbits  that  more  closely  follow  the  definitive  solutions 
based  on  additional  SLR  data.  However,  results  vary  depending  on  the  orbit  of  the  satellite  and  the  time 
past  the  initial  fit  span. 

Ultimately,  these  results  provide  a  relative  benchmark  of  the  geopotential  and  ocean  tide  models  for  the 
predictive  performance  of  precision  orbit  determination  solutions.  This  is  important  for  satellite  mission 
planners  who  often  require  accurate  predicted  orbits.  Future  work  includes  examining  other  models,  in¬ 
corporating  more  geodetic  satellites  and  examining  the  time  variability  of  these  results.  Additionally,  the 
comparison  metrics  in  this  analysis  provide  a  level  of  confidence  in  a  model’s  ability  to  capture  real  world 
dynamics  by  examining  the  consistency  between  a  predicted  orbit  and  a  fit  orbit  using  the  same  force  models. 
Future  research  can  expand  on  this  to  include  a  comparison  of  the  fit  orbit  to  an  ILRS  reference  orbit  as  a 
measure  of  accuracy. 


9  of  12 


American  Institute  of  Aeronautics  and  Astronautics 


VII.  Appendix 


Orbit  comparison  results  for  all  geopotential  and  ocean  tide  model  combinations  are  given.  Results  are 
generated  by  comparing  the  predicted  orbit  solution  to  the  fit  orbit  solution  for  the  desired  interval,  as  shown 
in  Figure  1.  Figure  6  shows  the  RMS  of  the  radial,  in-track,  and  cross-track  position  differences,  along  with 
the  RSS  position  difference  for  LAGEOS-1  over  each  five  day  interval  for  the  full  set  of  geopotential  and 
ocean  tide  model  combinations  examined  in  this  paper.  Figures  7,  8,  and  9  provide  the  same  information 
for  LAGEOS-2,  STARLETTE,  and  STELLA,  respectively. 

Figure  6:  LAGEOS-1  Difference  Between  Predicted  and  Fit  Orbits  -  Results  Summary 


FES2004  Position  Difference 

GOT48  Position  Difference 

TOPEX3  Position  Difference 

TOPEX4  Position  Difference 

Geopotential  Model 

Radial 

In-Track  Cross-Track 

RSS 

Radial 

In-Track 

Cross-Track 

RSS 

Radial 

In-Track  Cross-Track  RSS 

Radial 

In-Track  Cross-Track  RSS 

EGM96 

0.052 

0.361 

0.348 

0.504 

0.047 

0.319 

0.426 

0.535 

0.055 

0.377 

0.344 

0.513 

0.055 

0.349 

0.346 

0.494 

n 

EGM2008 

0.056 

0.366 

0.292 

0.471 

0.050 

0.323 

0.376 

0.499 

0.059 

0.382 

0.283 

0.479 

0.058 

0.354 

0.287 

0.459 

GIF48 

0.055 

0.359 

0.283 

0.460 

0.049 

0.317 

0.368 

0.488 
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0.375 
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0.467 

0.057 

0.346 

0.277 

0.447 

1 
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0.351 
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EIGEN6C3 

0.052 

0.364 

0.322 

0.489 

0.047 

0.322 

0.404 

0.519 

0.055 

0.380 

0.317 

0.498 

0.055 

0.352 

0.320 

0.479 

EGM96 

0.049 

1.556 

0.332 

1.592 

0.029 

1.409 

0.540 

1.509 

0.053 

1.604 

0.310 

1.634 

0.052 

1.505 

0.313 

1.538 

« 

EGM2008 

0.056 

1.572 

0.395 

1.622 

0.035 

1.424 

0.544 

1.525 

0.059 

1.620 

0.369 

1.663 

0.059 

1.520 

0.368 

1.566 

GIF48 

0.051 

1.555 

0.443 

1.618 

0.034 

1.407 

0.569 

1.518 

0.055 

1.603 

0.418 

1.657 

0.054 

1.503 

0.416 

1.561 

1 

GGMS05S 

0.051 

1.575 

0.341 

1.612 

0.031 

1.426 

0.526 

1.520 

0.055 

1.622 

0.315 

1.654 

0.054 

1.523 

0.316 

1.556 

EIGEN6C3 

0.050 

1.579 

0.334 

1.615 

0.030 

1.431 

0.526 

1.525 

0.053 

1.627 

0.308 

1.657 

0.053 

1.527 

0.310 

1.559 

EGM96 

0.084 

3.665 

0.577 

3.711 

0.042 

3.281 

0.874 

3.396 

0.084 

3.729 

0.540 

3.769 

0.085 

3.596 

0.539 

3.637 

To 

EGM2008 

0.090 

3.687 

0.757 

3.765 

0.050 

3.300 

0.937 

3.431 

0.091 

3.749 

0.722 

3.819 

0.091 

3.612 

0.718 

3.683 

GIF48 

0.081 

3.648 

0.838 

3.744 

0.041 

3.263 

0.986 

3.409 

0.082 

3.711 

0.805 

3.798 

0.083 

3.575 

0.801 

3.664 

1 

GGMS05S 

0.086 

3.699 

0.647 

3.756 

0.044 

3.312 

0.887 

3.429 

0.086 

3.761 

0.610 

3.811 

0.087 

3.628 

0.608 

3.680 

EIGEN6C3 

0.082 

3.709 

0.625 

3.762 

0.041 

3.323 

0.879 

3.437 

0.082 

3.771 

0.588 

3.817 

0.083 

3.635 

0.586 

3.682 

EGM96 

0.097 

5.557 

0.899 

5.630 

0.040 

4.761 

1.266 

4.926 

0.097 

5.701 

0.851 

5.765 

0.099 

5.489 

0.851 

5.556 

i 

EGM2008 

0.106 

5.559 

1.231 

5.694 

0.052 

4.760 

1.435 

4.972 

0.107 

5.706 

1.192 

5.830 

0.109 

5.492 

1.187 

5.619 

GIF48 

0.094 

5.493 

1.355 

5.658 

0.042 

4.699 

1.519 

4.938 

0.096 

5.640 

1.319 

5.793 

0.098 

5.426 

1.313 

5.584 

1 

GGMS05S 

0.100 

5.594 

1.046 

5.691 

0.044 

4.794 

1.326 

4.974 

0.101 

5.750 

1.002 

5.838 

0.103 

5.528 

1.000 

5.619 

EIGEN6C3 

0.095 

5.613 

1.005 

5.703 

0.040 

4.817 

1.305 

4.990 

0.096 

5.760 

0.960 

5.840 

0.098 

5.548 

0.957 

5.631 

EGM96 

0.072 

7.027 

1.014 

7.100 

0.018 

5.730 

1.449 

5.910 

0.060 

7.242 

0.963 

7.306 

0.067 

6.932 

0.954 

6.998 

n 

EGM2008 

0.075 

7.025 

1.445 

7.172 

0.010 

5.727 

1.651 

5.960 

0.066 

7.238 

1.408 

7.374 

0.071 

6.924 

1.391 

7.063 

GIF48 

0.059 

6.928 

1.599 

7.110 

0.009 

5.629 

1.753 

5.896 

0.051 

7.141 

1.566 

7.310 

0.056 

6.826 

1.548 

7.000 

1 

GGMS05S 

0.072 

7.088 

1.212 

7.191 

0.014 

5.789 

1.522 

5.986 

0.062 

7.302 

1.168 

7.396 

0.067 

6.989 

1.154 

7.084 

EIGEN6C3 

0.067 

7.120 

1.158 

7.214 

0.013 

5.821 

1.498 

6.011 

0.056 

7.333 

1.113 

7.417 

0.062 

7.018 

1.100 

7.104 

Figure  7:  LAGEOS-2  Difference  Between  Predicted  and  Fit  Orbits  -  Results  Summary 
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Figure  8:  STARLETTE  Difference  Between  Predicted  and  Fit  Orbits  -  Results  Summary 
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Figure  9:  STELLA  Difference  Between  Predicted  and  Fit  Orbits  -  Results  Summary 
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